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ABSTRACT 

We present results for the spectral distortions of the Cosmic Microwave Background 
(CMB) arising due to bound-bound transitions during the epoch of cosmological hy- 
drogen recombination at frequencies dow n to v ~ 100 MHz. We exten d our previous 
treatment of the recombination problem ( Rubi fkT-Martfn et al.ll2006|) now including 
the main collisional processes and following the evolution of all the hydrogen angular 
momentum sub-states for up to 100 shells. We show that, due to the low baryon den- 
sity of the Universe, even within the highest considered shell full statistical equilibrium 
(SE) is not reached and that at low frequencies the recombination spectrum is signifi- 
cantly different when assuming full SE for n > 2. We also directly compare our results 
for the ionization history to the output of the Recfast code, showing that especially 
at low redshifts rather big differences arise. In the vicinity of the Thomson visibility 
function the electron fraction differs by roughly —0.6% which affects the temperature 
and polarization power spectra by <J 1%. Furthermore we shortly discuss the influence 
of free-free absorption and line broadening due to electron scattering on the bound- 
bound recombination spectrum and the generation of CMB angular fluctuations due 
to scattering of photons within the high shells. 

Key words: atomic processes - cosmic microwave background - spectral distortions 
- cosmology: theory - early Universe 



1 INTRODUCTION 

The study of the Cosmic Microwave Background (CMB) 
provides one of the most powerful tools to test cosmologi- 
cal models and to determine the parameters describing the 
Universe. Nowadays, most of the CMB experiments are fo- 
cused on the study of the temperature and polarization an- 
gular fluctuations, which have been measu red recently with 
very high precision b y the Wmap satellite llPaee et alJl200(i 
Hins haw et all2 006*l. permitting to achieve accuracies in the 
determination of some of the key cosmological parameters 
on the level of a few percent. 

The CMB photons which are detected today were es- 
sentially last scattered towards us during the epoch of cos- 
mological recombination, where the temperature of the Uni- 
verse had beco me sufficiently low to p e rmit the forma tion of 
neutral atoms (IZeldovich et alJll968t iPeebled Il96sh . Dur- 
ing the process of cosmological recombination of the hy- 
drogen and helium atoms there was a net generation of 
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photons, which introduces distortions to the CMB black- 
body spectrum. Except for the far Wien tail of the CMB, 
these distortions are expected to be very small: the dis- 
tortions due to bound-bound transitions reach the level of 
A7/7 ~ 1CP 7 at frequency v ~ 1 GHz and are smaller 
than AI /I ~ 10~ 8 in the range 10 GHz < v < 500 GHz 
(e.g. see iRubino-Martm et alJ I2006L RMCS06 hereafter). 
This is mainly due to the extremely large entropy of the 
Universe (there are ~ 1.6 x 10 9 photons per baryon), 
since in the presence of undistorted blackbody radiation 
there cannot be much more than one photon per transition 
iLiubarskii fc Sunvaevlll983l) . The free-bound emission in- 
creases the spectral distortions by additional ~ 30—90% and 
a total ~ 5 photons per hydrogen atom are released due to 
hydrogen recombination in the undistorted ambien t black- 
body photon field of the CMB JChluba fc SunvaevlEx)6ah . 
However, a measurement of these tiny deviations from a pure 
blackbody would provide another independent way to deter- 
mine some of the key cosmological parameters, such as the 
baryon density or the total matter content. In addition they 
in principle are a byproduct of the accurate treatment of 
the recombination history and in particular the time depen- 
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dence of the free electron fraction in the Universe, which 
today is so important for the theoretical prediction of the 
CMB temperature and polarization power spectra. 

The numerical solution of the hydrogen recombina- 
tion history and the associated spectral distortions of the 
CMB requires the integration of a stiff system of cou- 
pled ordinary differential equations, describing the evolu- 
tion of the populations of the different hydrogen levels, 
with extremely high accuracy, simply because the ampli- 
tude of the distortions is so small. Several authors have com- 
puted the distortions due to bound-bound transitions for a 
mult i- level hyd rogen atom with different assumptions and 
simplifications toubrovicrJll97St Bcigman & Sunvacv 1978; 
Rvbicki fc dell'Antoniolll993l: iDubrovich fc Stolvarov 199E : 
Dubrovich fc Shakhvorostovj2004lKholuDenko et alJl200E : 
Wong et alJl2006l RMCS06). Among these the most impor- 
tant simplification is to assume full statistical equilibrium 1 
(SE) within a given shell for n > 2 (for a more detailed 
comparison of the different approached see RMCS06 and 
references therein). With this assumption the recombination 
problem becomes significantly simpler, since for each addi- 
tional shell only one mo re equation has to be treated, and 
iKhohioenko et all (j2005l> already presented results for the 
hydrogen recombination spectrum due bound-bound transi- 
tions in the frequency range 1 GHz ^ v 1000 GHz includ- 
ing up to 160 shells within this approximation. 

In our previous work (RMCS06) the spectral distortions 
of the CMB spectrum arising due to bound-bound tran- 
sitions and the 2s two-photon decay during the epoch of 
cosmological hydrogen recombination were obtained in the 
frequency range 1 GHz v ^ 3500 GHz. There the complete 
set of ordinary differential equations was solved for the cos- 
mological hydrogen recombination problem including up to 
30 shells, with all the energetically degenerate angular mo- 
mentum sub-levels treated separately. It was shown that as- 
suming SE for shells n > 2 leads to significant differences in 
the cosmological hydrogen recombination spectrum for the 
Balmer, Paschen and Brackett series. Furthermore, it was 
argued that collisions probably can only become important 
for shells well above n ~ 25 — 30, so that assuming SE for 
the lower shells is not justified, but a detailed examination 
of this problem has not been performed until now. 

In this paper, we present detailed results for the hydro- 
gen recombination spectrum due to bound-bound transition 
at low frequencies down to v ~ 100 MHz (see Fig. |2J. To 
this end, we have solved the complete set of equations for 
a multi-level hydrogen atom with 100 shells, i.e. we are fol- 
lowing the populations of 5050 separate hydrogen states. 
Furthermore, we include the main collisional processes into 
our calculations (for more details see Sect. 12. Q and discuss 
deviations from full SE within the high shells. We also exam- 
ine the influence of the free-free process (Sect. 14.11 at low 
frequencies, the line broadening due to electron scattering 
(Sect. R~2"jl and the generation of CMB angular fluctuations 
by scattering of photons within the high shells (Sect. l4~3t . 
For the treatment of the multi-level hydrogen atom we ba- 



1 i.e. the population of a given level (n, I) is determined by N n i = 
(21 + l)N n /n 2 , where N n is the total population of the shell with 
principle quantum number n. 



sically follow the same procedure as described in RMCS06, 
with some modification which are explained in Sect. 2. 

Several authors have worked on percent level corrections 
to the ionization history arisi ng from different previously 
neglected physical process es (|r2ubroyicIi_fc_Gj^£h^vl 2005; 
Chluba fc Sunvaevl l2006bl iKholupenko fc Ivanchikl bood: 
Novosvadl^n2006T7 . It was pointed out that these correc- 
tions lead to differences in the theoretical predictions for the 
CMB temperature and polarization power spectra on a sim- 
ilar level and therefore should be taken into a ccount when 
analyz ing future high precision data. Recently, iLewis et alJ 
(2006) have made some first steps towards quantifying the 
possible impact of percent level correction to the ionization 
history on the estimation of cosmological parameters, again 
showing that in the era of precision cosmology more efforts 
should be made to improve the calculation of the recombi- 
nation history. 

In RMCS06 it was also pointed out that there are cor- 
rections to the ionization history on the level of percent aris- 
ing due to the detailed treatment of the angular momentum 
sub-states of the hydrogen atom. With computations includ- 
ing up to 100 shells, one reaches the point to directly com- 
pare the obtained ionization history with the output of the 
Recfast code, which was developed to repres ents the com- 
putat ions based on a 300-level hydrogen atom iSeaeer et alJ 
1999), but under the assumption of SE for n > 2. We shall 
discuss our results for the obtained electron recombination 
history as compared to the Recfast output in Sect. 13.31 and 
show that especially at low redshifts rather big differences 
arise. However, still more work has to be done to reach per- 
cent level accuracy for the ionization history. 

For all the computation presented in this work the fol- 
lowing values of the cosmological parameters were adopted: 
fib = 0.0444, Qtot = 1, fim = 0.2678, fi A = 0.7322, 
Y p = 0.24, h = 0.71 and T = 2.725 K. 



2 COMPUTATIONAL ISSUES 

In general we use the same approach and definitions as in 
RMCS06. However, in order to solve the full problem for 
w max ^ 30 we had to improve the computational treat- 
ment significantly. Therefore we modified one of our pre- 
vious codes to include the most important collisional pro- 
cesses and made some more simplifying assumptions, which 
we shall briefly discuss in this Section. 

2.1 Inclusion of collisions 

It was pointed out in RMCS06 that the assumption of SE 
is equivalent to instantaneous redistribution of the electrons 
within a given shell. Especially for hydrogen levels (n, I) with 
large angular momentum quantum number (0 < K n - 1) 
this assumption changes the populations strongly, because 
for high shells due to the ^-dependence of the recombination 
Gaunt-factor direct recombinations to these levels become 
very unlikely. This point is illustrated in Fig. Q where we 
show the /-dependence of the recombination coefficients a n i 
for different shells 2 . One can clearly see, that for large n most 

2 Here o„j is defined such that N e iV p a ni is the change of the 
population of level (n, I) due to direct recombinations per sec- 
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of the electrons recombine to states with /// max ~ 0.1 — 0.2. 
Therefore in the full problem electrons mainly reach levels 
with large n and / connecting to other levels via radiative 
and collisional transitions, where the former are very ineffi- 
cient, since they proceeds via many intermediate states, and 
latter, due to the large entropy of the Universe, are rather 
slow. However, due to the Al = ±1 restriction for radiative 
dipole transitions, levels with large n and / quantum num- 
bers will depopulate rather slowly via low frequency tran- 
sitions which leads to a delicate interplay between all the 
rates connecting to the neighboring levels. 

In Fig. Q for n = 100 we also give the recombination 
coefficient within the Kramers' approximation. Here the de- 
pendence on the recombination Gaunt-factor is neglected 
such that the /-dependence of the recombination coefficient 
only enters due to the fact that one single electron can 
find (21 + 1) sub-levels to be captured to. Hence, within 
the Kramers' approximation electrons preferentially recom- 
bine to states with large I. On the other hand, it is clear that 
within the Kramers' approximation the photoionization pro- 
cess is also more effective for high /-states, which implies that 
after any disturbance for these levels full equilibrium with 
the continuum is reached on a much shorter time-scale. Since 
here we are looking for tiny deviations from full equilibrium 
therefore it is important to include the Gaunt-factors into 
the computations and to account for collisions which help 
changing the populations of high /-states. 

From Fig. Q one can also see that explicitly neglect- 
ing induced recombinations, but including the Gaunt-factor, 
significantly alters the /-dependence of recombination coef- 
ficient for large n. In addition the total recombination rates 
for high shells strongly decreases for this case. Below, we 
now outline which collisional processes were accounted for 
in the present paper. 



2.1.1 I -changing collisions 

One of the most important differences with respect to our 
previous computations (RMCS06) is that we now also in- 
clude collisions of the form (n,l) — + (n, I ± 1). These are 
important in order to restore and maintain full statisti- 
cal equilib rium (SE) within a given shell and from simple 
estimates ijPengellv fc Seatonlll964|) for the conditions like 
in our Universe one expects that within the redshift range 
500 ^ z ^ 2500 these start contributing for shells with prin- 
cipal quantum numbers n ;> 25 — 30. The corresponding 
coupling term for the population N n i of level (n, /) to the 
populations N n p of level (n, /') with /' = J± 1 due to the 
collision with a particle of species i is then given by 



dN nl 



dt 



(21 + 1) N nl , 
(21' + 1) N nl 



- 1 



(1) 



We calculate the collision rates C n i^ n y according to 
iBrocklehurstl lll97llh 

Since the energy of the bound electron is not changing, 
in contrast to n-changing collisions the mass of the projectile 
is more important than its velocity. Therefore the /-changing 
collision rates due to proton impact are usually larger than 



ond, where N e and N p are the free electron and proton number 
densities, respectively. 




Figure 1. /-dependence of the recombination coefficient, a n [, at 
z = 1300 for different shells. The curves have been re-scaled by 
the total recombination coefficient, o n ,tot = S; Q n!i an d multi- 
plied by /max = n — 1 such that the 'integral' over § = ///max 
becomes unity. Also the results obtained within the Kramers' ap- 
proximation, i.e. = const X [21 + 1], and without the inclusion 
of stimulated recombination for n = 100 are presented. 



for electron impact, but we included both into our calcu- 
lations. Also collisions with Hell ions were accounted for, 
but due to their low number density and the history of He II 
recombination they only lead to a < 10% contribution to 
the total collision rate at redshifts z > 1500. Since the lines 
due to hydrogen recombination typically appear at redshifts 
z < 1500, Hen collisions basically do not affect the final 
results for the bound-bound spectral distortion. 

In order to speed up our computation we tabulate the 
/-changing collision rates before starting to solve the system 
of coupled ordinary differential equations, assuming that a 
sufficient approximation to the electron temperature T e and 
electron, proton and Hell numb er densities is given by the 
results from the Recfast code iSeaeer et al]ll999l) . Since 
in general the approximations for the collisional rates are 
only accurate on the level of ~ 10% — 50% this should be a 
reasonable approach. 

One should mention that in principle the Al = ±1 re- 
striction does not apply for /-changing collisions and cou- 
pling terms with |A/| > 1 could be important, but smaller 
due to the decrease of the cross section with growing Al. 
These should help to bring the populations for a given shell 
faster into full SE than in the calculations presented here. 
However, we neglected those terms for the moment. 



2.1.2 n-changing collisions 

Energy or n-changing collisions of the type (n, /) — * (n , Z±l) 
with n^ti' are usually less important th an /-changing colli- 
sions (e.g. see comments in Sect. 3.2.4 of iHummer fc Storevl 
Il987l) . But since here we are discussing highly excited elec- 
tron states in the hydrogen atom we decided to include col- 
lisional excitation and de-excitation as well and check their 
influence on the recombination spectrum. However, we are 
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not including collisional ionization or three-body recombi- 
nation, since the relevance of these processes is expected to 
be even smaller. 

As mentioned above, for energy-changing collisions the 
velocity of the projectile is important and hence the n- 
changing collision rate is dominated by the contribution 
from electron impact. We shall only consider collisional exci- 
tation and de-excitation by electrons, but neglect the corre- 
sponding processes induced by protons, Hell ions or neutral 
hydrogen atoms. 

For collisional processes a vast amount of literature can 
be found, all based on different assumptions and approxi- 
mations. He re we decided to use the simple approximations 
as given by Ivan Regemorterl B1962I) . These likely are only 
accurate to within a factor of ~ 2 — 10 probably even up 
to ~ 100, but they are very fast for numerical evaluations 
and widely used in computations of stellar atmospheres. The 
corresponding coupling term for the population N n i of level 
(n, I) to the populations N n ni of level (n', I') due to the col- 
lision with an electron is then given by 



dN nl 



dt 



(21 + 1) N n , v 
(2Z' + 1) N nl 



A E 



(2) 

where AE is the energy difference of the considered levels. 
We calculate the collisional de-excitation rate C7 n !-m'!' f° r 
electron impact accor ding to formula (18) in the work of 
Ivan Regemorterl jl962l) . Only radiatively allowed transitions 
are taken into account here. Note that since collisions are 
controlled by electrons, the electron temperature appears in 
the exponential factor of 



2.2 Additional comments 

The computations of the photoionization and photorecom- 
bination rates involve 1-dimensional integrals over the am- 
bient photon field. These have to be evaluated many times 
in order to solve the problem. Since we assume that the am- 
bient photon field is a pure blackbody, the only coupling 
to the current solution arises due to the electron temper- 
ature. For up to 50 shells we have checked that instead of 
evaluating both the photoionization and photorecombina- 
tion rates at every time it is sufficient to use the photore- 
combination rate only and apply the detailed balance re- 
lation. Since this approximation is possible one can again 
tabulate the photorecombination rates using the results for 
the electron temperature from the Recfast code and then 
interpolate. In this way one speeds up the computation by a 
large factor. For 50 shells we have compared the results ob- 
tained with and without this simplification and found prac- 
tically no difference in the recombination spectrum. This is 
because at z ~ 1100 — 1400, i.e. where the lines form, the 
difference in the electron and photon temperature is very 
small (AT/T ~ 10~ 6 - 10~ 5 ). One should also mention that 
we evaluated the photoio nization cross section and radia- 
tive rates as described in IStorev ^ _lTrmimg^j[lggl 



4 

,yr 



for n > 50 the expressions given bv lKarzas fc Latterl ( 1961 
were numerically unstable. 

In addition it is important to use both relative and ab- 
solute error control. The solvers D02NBF and D02NGF from 



n = 100 

Spill 

: 180, n m =2 



V [GHz] 

Figure 2. The hydrogen bound-bound emission spectrum at low 
frequencies. In the full computation (ra max = fisplit = 100) the 
evolution of all the angular momentum sub-states was taken into 
account separately including collisions. For comparison we also 



show the results obtained with n„ 



180 and n 



split 



2. The 



contribution due to the 2s two-photon decay is negligible within 
the considered range of frequencies and was therefore omitted. 



the NAG-Library 3 provide this possibility. We usually re- 
quested e ~ 10~ 8 relative and e ~ 10~ 22 — 10~ 23 or less 
absolute error. Furthermore at low redshifts we had to limit 
the maximal step-size to h ~ 0.1 — 0.5. We also ran com- 
putations with increased accuracies but found no significant 
differences. With these settings one computation of the re- 
combination spectrum for n max = n sp iit = 100 took about 6 
days on a ~ 3 GHz single processor machine. 



3 RESULTS AND DISCUSSION 

3.1 The hydrogen bound-bound emission 
spectrum at low frequencies 

In Figure|5|we present the hydrogen recombination spectrum 
for bound-bound transitions at low frequencies. In our full 
computation (n max = n sp ]i t = 100) the evolution of all the 
angular momentum sub-states was taken into account sepa- 
rately and I- and n-changing collisions were included. In this 
calculation the evolution of 5050 separate hydrogen states 
in the redshift range of 200 ^ z ^ 2100 was followed. Due 
to computational limitations this was the largest calculation 
we were able to perform within reasonable time. For com- 
parison we also show the results obtained with n max = 180 
under the simplifying assumption of SE for shell with n > 2, 
i.e. n sp Ht = 2. 

For the full calculation (n max = 100 and n sp iit = 100) 
within the frequency range 1 GHz ^ v ^ 100 GHz one can 
clearly see the low contrast features, which were also dis- 
cussed in RMCS06. As expected, here due to the inclusion 
of more shells the level of the distortions is higher than in 
RMCS06, but at v < 2 GHz they do not reach the level ob- 
tained in the case n max = 180 and n sp nt = 2. At 1 GHz these 

3 see |http://www.nag.co.uk7| 
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Figure 3. The hydrogen bound-bound emission spectrum at low 
frequencies for different n max and n sp i; t = »i max . Also I- and n- 
changing collisions were included. 
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Figure 4. Effect of collisions on the hydrogen bound-bound emis- 
sion spectrum at low frequencies. Free-free absorption was not 
taken into account. 



two solutions differ by a factor of ~ 1.4 (see Section 13.1.21 
for more discussion). Below ~ 1 GHz the spectral distor- 
tions become practically featureless. In the frequency range 
1 GHz <J v <C 10 GHz the slope of the spectrum is /3 ~ 0.46 
and therefore slightl y steeper than the one obtained in the 
calculations of lKholupenko et al.1 <l2005h and our n max = 180 
and n sp iit = 2 case, i.e. /3 ~ 0.35. 



3.1.1 Convergence with n max 

To estimate the convergence of the results for the hydrogen 
recombination spectrum at low frequencies we computed the 
distortions for different n max (see Fig. [3J|. At v ~ 1 GHz for 
n max = 100 the spectral distortion is a factor of ~ 2.6 larger 
than for n max = 30. One can see that at v ^ 1 GHz the dis- 
tortion is practically not changing anymore when including 
more than ~ 80 shells. Within our assumption the results 



should be converged on a level of ~ 10 — 20% at frequencies 
0.1 GHz < v < 1 GHz and to better than 1% for v £ 1 GHz. 
We also checked the convergence for up to n max = 180 but 
with n S piit = 2 and reached to the same conclusion. 



3.1.2 The influence of collisions 

Here we wish to understand how the inclusion of collisions 
affects the hydrogen recombination spectrum. In Figurc^Jwc 
compare the results for the spectral distortions down to very 
low frequencies obtained for n max = 100 in three different 
calculations: case A: n sp i it = 2; case B: n sp ut = n max , no 



collisions; case C: n. 



put = n-n 



with collisions. 



First one can clearly see the difference in the slope 
of the distortions in the range 0.1 GHz — 10 GHz for the 
cases B and C in comparison with case A. At frequencies 
below v ~ 2 GHz the solution for case A lies well above 
the other two solutions. The solution for case B cuts off at 
v ~ 200 MHz, while the ones for the cases A and C reach 
down to v ~ 45 MHz. This shows that without collisions 
transitions with An < n at large n are disfavored and elec- 
trons reaching states 1 < n to a large fraction depopulate 
via transitions with An ~ n, which increases the emission 
at higher frequencies (see also RMCS06). In addition one 
can conclude that for case A and C the cascade of electrons 
within the high /-states starts to become similarly impor- 
tant. However, since the populations of the very low shells 
(n ^ 20 — 30) are not significantly altered by collision, at 
low frequencies the distortions for case C cannot reach the 
same level as in the case A, simply because the net radiative 
rates and population of the high shells critically depend on 
the cascade to lower shells. 

We have checked, that increasing the rates for l- 
changing collisions by a large factor (~ 10 4 — 10 6 ) we can 
recover the solution like for n sp ij t = 2. Furthermore we 
have performed additional calculations with n max > 100 
and n sp iit = 100 in order to check how this affects the re- 
combination spectrum. We found that it only influences the 
distortions close to the lowest frequencies and still there the 
level of the case A is not reached. Therefore we conclude 
that the difference in the slopes of the recombination spec- 
trum for case A and C at low frequencies will remain, even 
for n max = n S pi it > 100. 

We also studied how important in particular n-changing 
collisions are by performing computations only including l- 
changing collisions. In this case the spectrum practically re- 
mained unchanged (differences < 1% at the very low fre- 
quencies). In addition we ran computations artificially in- 
creasing the rate of n-changing collisions by a factor of 100. 
Here we found that the recombination emission at low fre- 
quencies was slightly reduced (at most by factor of ~ 2 for 
v < 100 MHz in this case). This is due to the fact that in this 
case the net rate connecting two levels with large n includes 
a sizeable non-radiative contribution. One would expect that 
this, in addition to free-free absorption, may lead to a re- 
duction of emission at very low frequencies, simply because 
n-changing collisions should dominate over radiative transi- 
tions for some large value of n, but it seems that this will 
only happen for n max > 100 and correspondingly at even 
lower frequencies than considered here. 
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Figure 5. Departures of the total populations from Saha: N n 
denotes the total population of shell n as obtain in our 100 shell 
computation with n max = n sp ij t including collisions. N n g a h a is 
the expected Saha-population relative to the continuum. 



3.2 Populations of the levels 

In this Section we discuss the structure of the hydrogen level 
populations in some more detail. At redshifts z 1500 the 
ground state population always is extremely far from Saha- 
equilibrium with the free electrons. We therefore omitted it 
in the following. 

Figure shows the departures of the total populations 
for given n from the expected Saha population relative to 
the continuum. At high redshifts [z > 1800), for all shells 
with n ^ 2 the departures from Saha are <J 10~ 7 — 10 -6 . 
At z ~ 1500 the departures reach a level of ~ 10 -5 — 10 -4 
and for z <J 1500 all the curves become flatter, because from 
there on recombination is controlled by the escape of pho- 
tons from the Lyman-a transition. At z ~ 800 the devia- 
tions from Saha exceed ~ 10% for all the considered shells 
and below that redshift they only change within a factor of 
~ 2 — 5 until redshift z = 200. In general one can also see 
that for shells with n < 10 the departures from Saha are 
typically larger than for higher n and that below z ~ 800 
the total populations of these levels are overestimated by 
many orders of magnitude when applying the Saha-relation 
with respect to the continuum. Note that during the whole 
epoch of recombination the total population of any shell is 
smaller than the expected Saha value. 

In addition one can look at the deviations of the popu- 
lations within a given shell from full statistical equilibrium 
(SE) over the Z sub-levels. In Figures |S| and |7| we show the 
non-equilibrium effects on the populations of the angular 
momentum sub-states at different redshifts and given n for 
our computations with n max = n sp iit = 100. In each panel 
the results obtained with and without the inclusion of Z- 
changing and n-collisions are presented. In general one can 
state that departures from full SE become smaller the higher 
the redshift and the larger the value of n is. For n = 30 the 
inclusion of collisions does not alter the results significantly 
at any redshift, while for n — 50 one can see that within 
the high Z-states the deviations from SE start to vanish. For 
n = 80 and n = 100 within the high /-states the deviations 
from full SE are very small down to redshifts z ~ 1100— 1200 
and even in the low Z-states departures from SE start to be 
wiped out. 

But still even for n max = 100 there is no shell beyond 
which full SE is reached for all the angular momentum sub- 
states. All the low Z-states are underpopulated as compared 
to full SE. This is likely due to the fact that for n > 40 these 
levels can all connect more or less directly with shells n < 30, 
which are always rather far from SE. Also the coupling to the 
continuum due to the Z-scaling of the photoionization Gaunt 
factor is much stronger than for large I. Even for n = 100 the 
/-changing collisions are not able to significantly modify the 
low-Z-populations set by the radiative rates. However, one 
still expects that increasing n max for some n > 100 also the 
low Z-states will reach full SE. Unfortunately due to com- 
putational constraints with n sp ut = n max we were not able 
to go beyond n max ~ 100. We conclude that for the cos- 
mological hydrogen recombination calculations Z-changing 
collisions become important for shells with n > 30 — 40, 
as already estimated in RMCS06, but n-changing collisions 
may be neglected for shells with n <, 100. 
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Figure 6. Non-equilibrium effects on the populations of the angular momentum sub-states for our computation with n ma x = ^split — 100 
at different redshifts and given n. The left column is for z = 900, the middle for z = 1000 and the right for z = 1100. The computations 
were performed for n ma x = 100 following the populations of all the angular momentum sub-states separately. In each panel the results 
obtained with and without the inclusion of /- and n-changing collisions are given. We present AN n i/N n i = [N n i — TV^Pj/iV®?, where the 
statistical equilibrium (SE) population is computed from the actual total population of the shell by N°? = [(2/ + l)/n 2 ]iVtot ■ 
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Figure 7. Non-equilibrium effects on the populations of the angular momentum sub-states for our computation with n max = "-split — 100 
at different redshifts and given n. The left column is for z = 1200, the middle for z = 1300 and the right for z = 1400. The computations 
were performed for n max = 100 following the populations of all the angular momentum sub-states separately. In each panel the results 
obtained with and without the inclusion of /- and n-changing collisions are given. We present AN n i/N n [ = [N n i — N~¥]/N??, where the 
statistical equilibrium (SE) population is computed from the actual total population of the shell by iV~p = [(2/ + l)/n 2 ]7Vtot- 
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Figure 8. Relative difference in the free electron fraction, 
AN e /N e , for different values of n max . The upper panel shows the 
comparison for our computations with n sp ]j t = n max with the 
results for n sp ii t = 2. In the middle panel we compare our com- 
putations for n sp Bt = 2 with the output of Recfast and in the 
lower panel we show the relative difference for our computations 
with n sp i; t = n max with respect to the Recfast output. 

3.3 Effects on the free electron fraction and 
comparison with Recfast 

In the upper panel of Fig. |H| we show the relative differ- 
ence in the free electron fraction for our computations with 
w S piit = Umax as compared to the results for n sp n t — 2. At 
low redshifts this difference exceeds the level of 5%, show- 
ing that when following the populations of all the angular 
momentum sub-states separately recombination is slightly 
slower. However, it seems that close to the maximum of the 



Thomson visibility function (zdec ~ 1090 and Az ~ 195) 
departures from full SE within the shells are not sufficiently 
large to modify the recombination history. At z ~ 500 the 
relative difference for the case n max = 100 is more than 2 
times bigger than for n max = 30. Also in the redshift range 
z > 700 the ratio of these two curves remains greater than 
~ 1.2. Including more than 100 shells into the calculation 
seems necessary in order to reach full convergence at all red- 
shifts. 

In the middle panel of Fig. |S]we present the direct com- 
parison of our results with the output of Recfast 4 for dif- 
ferent cases with n sp iit = 2 up to n max = 180. During the 
peak of the visibility function the difference is close to —0.7% 
and only reaches ~ 1% and ~ 3% for z ~ 960 and z ~ 780, 
respectively. This comparison shows that representing the 
full multi-level calculation with percent-level accuracy in the 
full range of redshifts probably requires some more general 
'fudge-function' instead of a single 'fudge-factor'. Also, the 
difference at z > 1000 is probably not related to the de- 
tailed treatment of the populations in the angular momen- 
tum sub-states but points towards the limitations of the 
effective 3-level approach. A detailed analysis will be left for 
some future work. 

Finally, in the lower panel of Fig. |H| we show the di- 
rect comparison of our results for different calculations with 
Umax = n sp iit with the output of Recfast. There is a dif- 
ference of roughly —0.6% in the electron fraction close to 
the peak of the visibility function. The results also sug- 
gest that at low redshifts recombination proceeds slightly 
slower than in the Recfast code and that one can still ex- 
pect more than a few percent corrections due to the detailed 
treatment of the populations in the angular momentum sub- 
states. However, one still has to push to a larger number of 
shells and very likely include more physics, such as details 
of helium recombination, collisional ionization and three- 
body recombination, collisional transitions with |A£| > 1 
or forbidden transitions. At low redshifts (z <; 300), also de- 
tails of the hydrogen chemistry may have to be reconsidered. 
In addition one should include a sim ultaneous treatment of 
other processes such as discussed bvlDubrovich fc Grachev 
|2005[) : IChluba fc Sunvael(l2006bl) : lKholupenko fc Ivanchik 
(2006); Novosvadlvi (2006) and likely more. Here especially 
the treatment of two-photon transitions for higher shells 
including stimulated emission or possibly other forbidden 
transitions may still lead to percent level differences in both 
the hydrogen recombination spectrum and the ionization 
history. 

As an example, the strongest distortions of the CMB 
spectrum arising from the epoch of hydrogen recombination 
are due to the Lyman-a transition and the 2s two-photon 
decay. Previously, the feedback of these excess photons onto 
the photoionization rates for the secon d shell (n — 2) has 
been considered bv ISeaeer et alJ fcOOOh . yielding no signifi- 
cant changes in th e ionization history. However, r ecently it 
has been shown by Kholupcn ko fc Ivan chik (2006) that this 
huge excess of photons in the Wien-tail of the CMB leads 
to an increase of the Is — > 2s two-photon absorption rate, 



4 In fact we used our own implementation of this code, 
which solves the system of differential equations as given in 
ISeager et al.lll999l) with accuracy better than 0.1%. 
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Figure 9. changes of the CMB temperature (TT) and polariza- 
tion {EE) power spectra for our full computation with n max = 
n split = 100 including collisions. 

which results in percent-level corrections to the ionization 
history during hydrogen recombination. This feedback can- 
cels a large part of the effect o f induced two-photon decay 
of the 2s level as discussed bv IChluba fc Sunvaevl (|2006b) 
which again shows that reaching percent-level accuracy for 
recombination is not very simple. 

We would like to point out that the simplification in 
connection with the computation of the photoionization and 
recombination rates, as described in Sect. 12.21 at low red- 
shifts may still affect the results presented in Fig. |H]on the 
percent-level. 

3.4 Changes in the CMB temperature and 
polarization power spectra 

In order to compute the changes of the CMB temperature 
and polarization power spectra due to the differences in 
the ionization history we modified t he Recfast routine in- 
side the latest version of Cmbeasy <IDoranll2005l) . Since we 
obtained the ionization history only in the redshift range 
200 ^ z < 2100 we had to merge the solution of Recfast 
with our solution. At z = 2100 this was not problematic 
since there both solutions practically coincide. At z ^ 200 
we simply rescaled the Recfast output by a constant fac- 
tor, such that the merged solution became continuous. We 
checked how this approach influences the obtained difference 
in the power spectra by changing the transitions redshift to 
2 = 300. However, this only affected the changes in the 
power spectra at I <J 50 — 100 on the level of percent and 
we expect that the results for I .> 100 — 200 are robust. 
Note that here / corresponds to the index of the spherical 
harmonic function Yi m rather than the angular momentum 
quantum number for a level inside the hydrogen atom. 

In Fig. |5]we give the difference in the temperature and 
polarization power spectra for our full computation with 
n max = n sp iit = 100. The difference reaches the level of 
0.4% for the TT and 0.8% for the EE power spectrum. We 
checked the convergence by using the results for different 
n max < 100 and found that the curves still changed at the 




V [GHz] 

Figure 10. Effect of free-free absorption on the hydrogen bound- 
bound emission spectrum for n max = 180 and n sp i it = 2. The 
spectrum in the range u > 100 MHz is not affected by the free- 
free process, whereas at v < 1 MHz all the distortions are erased. 

percent-level when going from n max = 80 to 100. This shows 
that still more shells should be included and computations 
have to be extended down to z = in order to obtain a 
definite answer. 



4 ADDITIONAL PHYSICAL PROCESSES 

4.1 Effects of free- free absorption 

The free-free process is increasingly important at low fre- 
quencies and should eventually be able to wipe out any 
spectral distortion below some frequency vr. Given the solu- 
tions for the recombination history and the emission result- 
ing from each bound-bound transition one can also take the 
free- free process into account. Since at redshifts z > 1200 
the difference in the electron and photon temperature is 
very small, for simplicity we assumed T c — T 7 = To(l + z) 
with todays CMB temperature T = 2.725 K. In Figure ITul 
we illustrate the effect of free-free absorption on the hydro- 
gen recombination spectrum for n max = 180 and n sp n t = 2. 
We chose this particular example, since we could obtain the 
distortions down to sufficiently low frequencies. Our calcula- 
tions show that the recombination spectrum in the frequency 
range v > 100 MHz is not affected by the free-free process. 
However, for He ill and probably Hell recombination the 
effect of free-free absorption will be stronger, since the tem- 
perature and density of the plasma is larger and there are 
more electrons per helium atom when hydrogen is still fully 
ionized. 

4.2 Effects of electron scattering 

To estimate the effects of electron scattering one can calcu- 
late the expected broadening of an initial ly narrow line for 
a give n recombinations history. Following IPozdniakov et alj 
( 1979) for large Thomson optical depth r = J N b (jt dl ^> 1 
and small Compton y-parameter y = J " 5 dr < 1 the 
broadening due to electron scattering is dominated by the 
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Figure 11. Line optical depth for some a-transitions within high 
shells. The results were obtained using the computations for 100 
shells. For comparison the H a -line (3 — ► 2) optical depth is given. 

Doppler effect and one can find: 4r| sc ~ 2VyTn2. Using 
the standard ionization history as obtained with the Rec- 
fast code one finds that for the spectral distortions due to 
hydrogen recombination the broadening is less than percent 
for emission occurring at redshift z < 1500. Since the feature 
we are considering here are typically much broader (> 20%) 
one may neglected the effect of electron scattering on the 
hydrogen recombination spectrum. However, for Hell and 
in particular for He ill recombination this may not be true 
anymore due to two reason: first He ill recombines at higher 
redshift (z ~ 6000-6500) than hydrogen (z ~ 1100-1500) 
and second because it occurred more rapidly. We shall not 
discuss these points any further here. 

4.3 Imprint of the recombinational line transitions 
on the CMB angular power spectrum 

Our detailed study of the hydrogen r ecombination permits 
us to evaluate the effects discussed on lRubino-Martfn et alJ 
(2005), but now for shells above n = 10. In that pa- 
per, these authors discussed the imprint of the recombina- 
tional lines on the angular power spectrum of the C MB, 
via resonant scattering in the lines (see also iBasu et alJ 
l2004l : lHernandez-Monteagudo fc Sunvaevl2005l . for more de- 
tails). Due to correlations with the existing radiation field, 
the net effect on the angular power spectrum, SCe, is pro- 
po rtional to the optical dept h in the line, tl (see Eq. 1 
in iRubino-Martm et alj |2005[). at the relevant redshift. As 
shown in iRubino^Martmei^al] j2005l) . in the optically thin 
limit the first order correction to the visibility function has 
two terms. The first one is connected with suppression of ex- 
isting anisotropics, and it is proportio nal to tl x Vt, where 
Vt is the Thomson visibility function iSunvaev fc Zeldovicbl 
ll970T) . The second term is associated to generation of new 
anisotropies, and it is proportional to tl x exp (— tt). Note 
that here I corresponds to the index of the spherical har- 
monic function Yj m rather than the angular momentum 
quantum number for a level inside the hydrogen atom. 

Since the bound-bound recombination lines are very 



narrow as compared to the duration of recombination, fix- 
ing the observing frequency one tunes to the signal arising 
during one particular epoch of recombination. Under the as- 
sumption that the maximal signal for all the bound-bound 
recombination lines arises at the same redshift this implies 
that going to lower frequencies the effect is dominated by 
the signal coming from the a-transitions of higher shells. 

In Fig. llll we present the optical depth to resonant scat- 
tering for several main transitions (An = 1) between high 
shells of the hydrogen atom and we compare these to the one 
corresponding to H a . Although the optical depths for n J> 20 
exceed the one for H Q at high redshifts, the expected signa- 
ture in the power spectrum will be small due to the presence 
of Thomson scattering. However, this process might become 
interesting in the future because at low frequencies the sig- 
nals of many transitions overlaps and one expects that for 
a given bandwidth of the experiment a superposition of the 
contributions with comparable amplitude might amplify the 
signal. This may increase the total change in the tempera- 
ture power spectrum notably, especially at frequencies be- 
low ~ 1 GHz. We have computed the effect on the power 
spectrum and found that for v ~ 50 — 100 MHz it is com- 
parable to the signal produced by the scattering within the 
H a -transition. Our computations also showed that a similar 
process due to the absorption and re-emission of photons in 
the continuum also leads to weak fluctuations on the level 
of 5C t l{l + 1)/2tt « 9 x 10" 4 uK 2 at 1.4 GHz. This signal 
also shows a characteristic oscillatory behaviour, and again 
it is maximum at small angular scales (£ ~ 900). 



5 CONCLUSION 

We presented results for the spectral distortions of the CMB 
arising during the epoch of cosmological hydrogen recom- 
bination due to bound-bound transition with observing fre- 
quencies down to v ~ 100 MHz and confirm that the relative 
spectral distortions, AI V /B V , are rapidly growing towards 
low frequencies. We show that even within the highest con- 
sidered shell full statistical equilibrium is not reached and 
that at low frequencies the recombination spectrum is sig- 
nificantly different when full SE for shells with n > 2 is 
assumed. We also conclude that in cosmological hydrogen 
recombination calculations /-changing collisions should be 
included for shells with n J> 30 — 40, but n-changing colli- 
sions can be neglected for shells with n < 100. 

In Sect. l3.3l we have presented our results for the ioniza- 
tion fraction as compared to the Recfast output. Including 
a treatment of the populations of all the individual angular 
momentum sub-states slows recombination down, leaving a 
> 5 — 10% larger residual electron fraction at low redshifts. 
Although the overall effect on the temperature and polar- 
ization power spectra due to the changes in the recombi- 
nation history is < 1% these additional electrons may be 
important in computat i ons of the first stars and structures 
jBarkana fc Loebll200ll iBromm et alJl2002ft . for which es- 
pecially the formation of molecular hydrogen plays a cru- 
cial role. The discussion has shown that representing the 
full multi-level calculation with percent level accuracy in 
the full range of redshifts probably requires some more gen- 
eral 'fudge-function' instead of a single 'fudge-factor'. Also, 
one still has to push to a larger number of shells and very 
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likely include more physics. For extensive calculations of the 
CMB temperature and polarization power spectra it will be- 
come important to develop a new simplified and fast scheme 
which accurately models the physics during the recombina- 
tion epoch. 
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